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electron  transfer  to  the  incoming  projectile  ion,  the  Anderson  correlation 
energy  U  is  included  in  the  analysis.  Calculations  of  the  probability  of 
negative-ion  formation  are  in  qualitative  agreement  with  experiments  on  the 
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I.  Introduction 


This  paper  addresses  recent  advances  which  we  have  made  in  the  theoretical 
description  of  charge  transfer  occurring  in  ion-surface  scattering,  where  one  or 
two  electrons  are  transferred  from  the  target  surface  to  the  incoming  projectile 
positive  Ion.  In  Section  II  we  discuss  the  transfer  of  one  electron  from  target 
surfaces  to  incoming  singly-charged  ions.  We  start  by  giving  a  phenomenological 
description  of  two  major  features  resulting  from  an  ion-neutralization 
experiment  in  which  the  neutralization  probability  P  versus  the  projectile 
collision  energy  E  is  detected:  (i)  exponential  decrease  of  P  with  the  increase 
of  E  in  the  high  collision  energy  region  and  (ii)  oscillatory  behaviour  of  P(E) 
in  the  low  collision  energy  region.  We  proceed  by  outlining  a  perturbati ve 
approach  to  the  chosen  ion-surface  interaction,  stating  the  rules  to  compute  at 
any  order  in  the  perturbative  expansion  the  neutralization  probability  and  by 
analyzing  the  conditions  under  which  the  expansion  converges  rapidly.  Finally, 
we  present  calculations  of  the  neutralizat ion  probability  P(E)  of  protons  from 
various  alkali-halide  surfaces.  For  these  systems,  the  fast  convergence 
conditions  are  met,  and  the  dependence  of  P  on  the  electronic  structure  of  the 
collision  partners  and  on  the  dynamics  of  the  process  are  critically  discussed 
at  the  end  of  Section  II. 

The  transfer  of  tw^  oi^ptrons  that  leads  to  the  negative-ion  formation  in 
ion-surface  scattering  is  discussed  in  Section  III.  Within  the  framework  of  the 
time-dependent  An dcrc :r model,  we  consider  the  two-electron  transfer 
process  in  ion-surface  scattering,  taking  into  the  account  the  Anderson 
correlation  energy  U  between  the  two  electrons  of  opposite  spin  in  the  same 
ionic  energy  level.  The  probability  of  negative-ion  formation  is  calculated  by 
means  of  the  time-evolution  operator.  The  results  are  in  qualitative  agreement 
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with  experiments  on  the  conversion  of  H  (D  )  to  H  (D  )  by  scattering  from  a 
W(  1 1 0 )  surface.  This  unified  point  of  view  is  further  extended  to  study  the 
temperature  effects  on  negative-ion  formation  from  surface  scattering.  The 
interplay  between  temperature  and  the  correlation  energy  U  is  discussed,  along 
with  new  ideas  for  improving  the  calculations  by  incorporating  the  image 
potential  into  the  ionic  energy  level  shift.  Finally,  Section  IV  is  the 
summary. 

II.  One-Electron  Transfer 

A  multiplicity  of  events  are  often  involved  in  ion  scattering  from  surfaces, 
where  the  neutralization  of  the  ion  is  one  of  them.  The  formation  of  neutral 
backscattered  particles  is  studied  by  monitoring  the  charge  state  distribution 
of  the  outgoing  particles  and  measuring  the  ratio  of  the  scattered-neutral  flux 
to  the  scattered-ion  flux.  When  more  information  is  required  on  the  precise 
state  of  the  neutralized  species,  then  the  photon  emission  from  them  is 
monitored.  The  mechanism  of  the  charge-transfer  process  must  be  understood  in 
order  to  explain  the  behavior  of  the  neutralization  probability  as  a  function  of 
the  collision  energy,  collision  geometry,  final  state  of  the  atom,  nature  of  the 
solid  target,  etc. 

A  typical  result  from  the  experiment  of  ion  neutralization  at  surfaces  shows 

the  following  features:1  The  neutralization  probability  P  as  a  function  of  the 

collision  energy  E  exhibits  a  maximum  at  some  value  E.  For  values  E  >  E,  P(E) 

manifests  a  rapid  decrease  to  zero,  and  for  E  <  E  oscillations  are  often 

displaced.  A  popular  phenomenological  description  of  the  process  is  summarized 
2-H 


below 


Let  us  assume  that  an  ion  A  is  scattered  from  a  surface,  where  the  atomic 
center  of  mass  follows  a  classical  trajectory  R(t).  The  total  Hamiltonian  of 
the  system  can  be  written  as 

H  «  HQ  +  HjCR(t)]  ,  (1) 

where  is  the  perturbation  created  by  the  incoming  ion.  If  R(t)  is  known,  H ^ 
becomes  a  function  of  time  t.  Let  Y<t)  be  the  rate  with  which  an  electron  is 
transferred  from  the  solid  to  the  ion  at  time  t.  Therefore,  setting  y  *>  1  -  P, 
we  can  write 


y(t) 

*  -Y(t)  y(t), 

(2) 

fl 

y(t) 

=  exp  |-l  dt  "Y ( t ) }  , 

(3) 

J0 

and  if  the  projectile  moves  perpendicular  to  the  surface  with  constant  velocity 
v  along  the  x-axis  we  have, 

2  r 

y(t-»®)  =  exp|-  -  |  dx  Y(x)}  ,  (U) 

J0 

where  x  ■  v|t|.  Assuming  also  that  Y(x)  «  YQe  ax,  we  get 


The  oscillations  are  seen  to  originate  from  quasi-resonant  electron  interchange 

5  6 

with  the  core  levels.  '  In  fact,  if  we  examine  a  two-(resonant)state  time- 


dependent  coupling,  the  Hamiltonian  matrix  can  be  written  as 


H 


V  ( t )  e 


V(t)e~i<Ht)  u> 


i  4>(  t ) 


(7) 


If  4>(  t )  is  assumed  to  be  independent  of  t  so  that  there  is  no  need  to 
considering  time-ordering  in  computing  the  evolution  operator,  and  if  the  phase 
$  of  the  interaction  is  taken  to  be  zero,  then  the  state  vector  of  the  system  is 


<K  t)> 


0>  cos 


dt'V(t’)]  -  | 1 >  i  sin 


J-C 


dt'V( t' )  ] 


(8) 


where  the  zero  of  the  energy  has  been  set  at  the  value  w,  and  where  |»p(-®)>  » 
1 0> .  The  probability  of  having  |il>(“’)>  -  jl>  is 


|<1  |  U  ( “ )  >  | 2 


sin2  [ 


dt’V(t')] 


or 


r 

sin2  I  dR  V ( R )  1 
•*0 


(9) 


which  is  an  oscillating  function  whose  amplitude  decreases  with  increasing 
collision  velocity  v. 

A  model  calculation  by  Bloss  and  Hone,^  which  includes  interactions  between 

broad-band  delocalized  electrons  and  localized  core  electrons,  exhibits  the 

characteristic  oscillations  (overimposed  on  the  exponential  decay).  The  Hamil- 
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tonian  used  is  an  Anderson-Newns  type  ’  to  which  is  added  a  core  level  |0> 

of  the  surface  atom  and  an  additional  state  |2>  on  the  ion: 

H  -  lc 
k 

+  I[Vk1(t)c*Cl  +  h.c.  ],  (10) 

k 


kCkCk  +  e0C0C0  +  e1C1C1  +  Z2 C2°2  *  ^V01(t)c0C1 


h.c.  ] 


where  k  labels  the  surface  delocalized  electrons.  The  following  items  are 
excluded  from  the  model:  (i)  electron-electron  interactions  (therefore  Auger 


1 


,  .j. 


neutralization  and  Auger  deexcitation);  (ii)  radiative  deexcitation;  (iii) 
negative-ion  formation  (therefore  Coulomb  correlation — Coulombic  repulsion  has 
been  taken  to  be  infinite);  (iv)  spin  dependency  [consequence  of  (iii)];  (v) 
image  forces  (therefore  the  atomic  energy  levels  have  been  taken  as  constants);  and 
(vi)  overlap  integrals  of  atomic  and  surface  states,  which  have  only  a  small 
contribution. 

il  6  Q 

Although  several  authors  ’  '  took  the  approach  of  solving  the  Heisenberg 
equations  of  motion  for  the  creation  and  annihilation  operators  c*  and  c, 
another  possible  approach  is  the  use  of  diagrammatic  techniques10  which  are 
widely  applied  to  time-dependent  problems.11  Due  to  its  simplicity,  the  former 
approach  is  more  commonly  used.  On  the  other  hand,  a  perturbative  treatment 
(with  diagrams)  allows  for  a  systematic  series  of  approximations  to  many-body 
problems,  and  we  shall  consider  this  below. 

We  consider  the  following  Hamiltonian: 

"u)  '  l  SVk  *  Wa  *  Hi(t)  • 

k 

where 

HI(t)  “  E  [Vk(t)cdck  +  (12) 

k 

Vk(t)  -  l/k  e”*!1!  ,  (13) 

and  X  is  a  parameter  measuring  the  duration  of  the  interaction.  A  single¬ 
band/single-atomic  level  interaction  is  considered,  and  we  want  to  explore  the 
consequence  of  a  perturbation  expansion  of  the  diagonal  elements  of  the  reduced 
density  matrix  of  the  neutralized  atom.  The  vacuum  state  is  the  state  in  which 
all  the  band  levels  |k>  are  filled  with  electrons  up  to  the  Fermi  energy  and  the 
ion  is  in  its  ground  state  so  that  the  state  |d>  is  unoccupied.  We  indicate  the 
state  of  the  system  in  which  one  electron  has  Jumped  from  |k>  into  jd>  by  | 1 d> 

-  cj(t)h^(t)|f  >,  where  1 1|/  >  is  the  (unknown)  ground  interacting  vacuum  state  in 


the  Heisenberg  picture  and  and  are  holt  ladder  operators.  The  charge- 
transfer  probability  is  therefore 


^oiy^d^iy 

<*01*0  > 


(ID) 


where  the  sum  over  k  is  nothing  more  than  a  partial  trace  over  the  state  of  the 

unobserved  surface.  In  the  interaction  picture  (in  which  operators  are  denoted 

oelow  by  a  tilde),  the  neutralization  probability  P  can  be  written  in  terms  of 

12 

the  (known)  non-interacting  ground  vacuum  state  |0>  as 


r 


yi  v 


<0 1 T | exp [ - 1 1  dt  e  hkc  V( t )  ]  1 1 0> 


<0 |T jexp[-i 


(15) 


dt  e"Vltlv(t')]!|C> 


where  V(t)  is  V ( t )  of  Eq .  (13)  in  the  interaction  picture. 

We  have  computed  the  neutralization  probability  P  up  to  third  order  in  the 

12 

perturbative  expansion  for  several  values  of  the  parameter  A  [Eq.  (13)].  This 

is  the  only  parameter  which  controls  the  dynamics  of  the  collision  process;  in 
particular,  it  controls  the  duration  of  the  interaction  and  hence  the  collision 
energy.  We  have  found  that  the  first-order  term  gives  the  main  contribution  to 
P  provided  (i)  the  duration  of  the  interaction  is  short  enough  or  (ii)  the 
atomic  discrete  level  is  not  embedded  in  the  valence  band  of  the  solid  where  the 
neutralizing  electron  originates.  The  maximum  exhibited  in  P(E)  in  typical 
experiments1  is  interpreted  to  result  from  the  interplay  between  two  factors: 
the  resonance  and  the  duration  of  the  interaction  factors.  As  the  energy 
increases,  the  former  makes  the  neutralization  easier  because  the  resonance 
condition  becomes  less  stringent,  but  the  latter  makes  the  neutralization  more 
difficult  because  the  electron  has  less  time  to  "jump"  into  the  ion.  Another 
factor  which  determines  the  magnitude  of  the  neutralization  probability  is  the 
relative  position  between  the  atomic  level  and  the  valence  band  of  the  solid. 
Excluding  the  case  in  which  the  atomic  level  is  embedded  in  the  band  (a 


8 


situation  in  which  the  first  few  terms  of  the  perturbation  expansion  are  not 
sufficient  to  give  a  converging  neutralization  probability),  the  closer  the 
discrete  level  is  to  the  continuum  of  the  band,  the  easier  the  charge  transfer 
occurs.  Finally,  a  wider  valence  band  always  enhances  the  neutralization 
probability. 

It  is  the  interplay  of  all  of  the  above  general  rules  that  determines 
quantitatively  the  relative  neutralization  probability  from  different  solid 
systems  into  different  atomic  states.  Therefore,  a  glance  at  the  electronic 
structure  of  the  atom-solid  pair  is  often  very  useful  in  determining  qualitative 
trends  in  the  neutralization  probabilities.  As  an  example,  we  display  in  Fig.  1 
the  electronic  structure  of  the  first  four  hydrogenic  shells  and  of  the  valence 

band  of  some  alkali-halide  solids.  Since  these  systems  have  already  been 

13 

discussed  at  length,  we  only  want  to  point  out  that,  for  example,  both  the 
width  of  the  bromide  valence  bands  and  their  proximity  to  the  hydrogen  ground 
state  would  (according  to  the  above  rules)  favor  the  trend  LiBr  >  NaBr  >  KBr 
for  all  values  of  the  interaction  times.  This  is  indeed  what  our  calculations 
have  shown. ^  The  calculation  of  the  neutralization  into  the  first-excited 
atomic  state  from  the  fluoride  systems  requires  a  more  exact  treatment.  At  very 
high  values  of  A,  where  the  resonance  factor  does  not  have  any  role,  the  valence 
band  width  determines  the  trend  LiF  >  NaF  >  KF ;  at  low  values  of  \  where  its 
increase  makes  the  neutralization  more  probable  due  to  the  ease  of  the  resonance 
condition,  the  proximity  of  the  n-2  level  to  the  upper  edge  of  the  valence  band 
determines  the  trend  NaF  >  KF  >  LiF  (see  Fig.  2).  The  absence  of  the  quasi¬ 
resonant  interaction  between  the  ionic  level  and  discrete  inner-core  surface 
states,  and  of  the  inclusion  of  many  terms  in  the  perturbative  expansion  for 
ionic  levels  embedded  into  the  solid  band,  precludes  the  observation  of  typical 
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oscillations  on  the  low-energy  side  of  the  neutralization  probability  plot.  The 
inclusion  of  long-range  image  forces,  and  therefore  of  a  more  realistic  time 
dependence  of  the  ionic  levels,  has  been  also  neglected.  We  shall  address  this 
aspect  in  the  analysis  of  the  two-electron  charge-transfer  process  at  the  end  of 
the  next  section. 

Ill .  Two-Electron  Transfer 

For  positive-ion  scattering  from  a  metal  surface,  there  is  considerable 

experimental  evidence  for  two-electron  transfer14  as  well.  In  experimental 

studies  of  two-electron  surface  scattering  processes,  for  example,  positive  ions 

such  as  H  ,  D  and  Na  are  scattered  from  a  W ( 1 10)  surface.  As  the  ions  leave 

the  surface,  a  certain  fraction  of  them  are  converted  into  negative  ions.  For 

the  purpose  of  increasing  this  fraction,  an  overlayer  of  cesium  is 

19  20 

deposited  on  the  metal  surface  to  reduce  the  work  function.  ’  The  formation 

of  negative  ions  from  surface  scattering  has  recently  attracted  more  interest 

21 

due  to  its  possible  use  in  plasma  diagonistics,  the  production  of  high 

22  2 3 
energetic  neutral  beams  and  heating  fusion  plasmas. 

Although  theoretical  investigations  have  been  made  successful  for  one- 

electron  transfer  processes,  ’  ^  ^  ’ 1  1  ^  based  on  the  time-dependent 

7  & 

Anderson-Newns  model,  the  theory  for  two-electron  transfer  processes  has 

27 

been  basically  phenomenological.  As  mentioned  before,  in  the  application  of 
that  model  to  the  process  of  ion  neutralization,  the  correlation  energy  due  to 
the  Coulombic  repulsion  between  the  two  electrons  of  opposite  spin  in  the  same 


ionic  energy  level  is  neglected  or  implicitly  assumed  to  be  infinite.  Since 
such  correlation  energy  has  been  shown  to  play  a  crucially  important  role  in 


localized  magnetization  around  impurities  in  transition  metals,  in  the  degree 
of  neutrality  of  hydrogen  atoms  chemisorbed  on  transition  metal  surfaces  and 


■  ■  •-  ■-  • 
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In  mixed  valences  occurring  in  rare-earth  compounds,  °  we  have  decided  to 

29 

consider  negative-ion  formation  by  taking  into  account  the  finite  Anderson 

7  8 

correlation  energy  U  in  the  time-dependent  Anderson-Newns  model. 

The  Hamiltonian  under  consideration  for  the  negative-ion  formation  is, 


H { t )  -  Hq  +  HI ( t )  , 


where 


r  t  r  + 

H  =  )  c  c,  c,  +}.  c  .c  .  c  .  +  U  n,  n  . 

0  ,  kkoko  ddodo  do  d,-o 

ko  o 


HI(t)  “  *  ^Vk(t)cdocko  +  h-C*^  ‘  (18) 

k,  o  "  ' 

V ^ ( t )  is  given  by  Eq.  (13).  and  n*o  =  CdoCdo"  Tbe  Probability  negative-ion 
formation  at  time  t,  P(t),  can  be  calculated  with  he  help  of  the  time-evolution 
operator,  T(t).  To  lowest  order,  P(t)  is  given  by 


dr  (  dr'  p(e)  p(e’)  <k+k| |T (2) (t) | I> |2  , 


whe.  Cp  and  tL  are  the  Fermi  energy  and  lower  conduction  band  edge  of  the 

metal,  respectively,  p(e)  is  the  density  of  the  conduction  band  and  is  assumed  to  be 

_ .. _ _  .1 _ _ »  .1 _ V,  J  6,9,25  ,L.  J  „  i  t  {  =  f  i.,^1  I  (S  -....-I  I  I,  1,  I  X 


constant  throughout  the  band, 


the  initial  and  final  states  |i>  and  |k+k^>, 


are  explained  in  Fig.  3,  and  f  y(t)  is  the  second-order  time-evolution  operator 
T(t)  in  the  interaction  picture.  Assuming  that  t  -  -*  when  the  incident  ions 
are  infinitely  far  away  from  the  surface,  we  have 


~  (  p  )  p  1  *  '  - 

r  ;(t)  -  (-1)  I  dt'  Hjit')  I  dt"  Hj ( t" )  , 

*  —  0D  *  —CD 


where 
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Figure  Captions 

1.  Electronic  structure  of  the  first  four  hydrogenlc  shells  and  of  the  valence 

band  of  some  alkall-hallde  solids.  Except  for  charge-transfer  from  the  LiF 

surface  into  the  hydrogen  ground  state,  the  relative  qualitative  behavior  of 

the  neutralization  probabilities  from  various  solids  into  various  atomic 

bound  states  can  be  predicted  by  comparing  the  relative  distance  between  the 

atomic  levels  and  the  valence  bands  with  the  relative  width  of  the  latter 

and  by  taking  into  account  that  the  frequency  of  the  n-th  shell  is 
2 

proportional  to  n  . 

2.  Neutralization  probability  P  into  the  first  hydrogenic  excited  state  (n-2) 
from  alkali-fluoride  surfaces  as  a  function  of  X.  Comparing  with  Fig.  1  , 
this  behavior  is  easily  understood:  both  the  resonance  and  bandwidth  factors 
determine  that  P(NaF)  >  P(KF).  At  low  values  of  X  where  the  resonance 
factor  is  important,  P(KF)  >  P ( LiF) ;  at  high  values  of  X  where  the  bandwidth 
factor  is  important,  P(LiF  >  P(NaF). 

3.  Schematic  diagram  of  electron  transfer  in  negative-ion  formation  at  T  -  OK. 

is  the  upper  and  eL  the  lower  edge  of  the  band,  (a)  Initial  state  |l>: 

valence  state  ed  empty,  band  filled  up  to  the  Fermi  level  ep.  (b) 

Intermediate  states  |k+>  and  |k}>:  states  corresponding  to  the  neutralized 

atom;  one  electron  transferred  from  the  e,  or  e,'  level  to  the  eJ  level.  The 

k  k  d 

arrows  denote  the  spin  states  of  the  -ileclrons,  ana  the  solid  and  hollow 
circles  represent  electrons  and  holes,  respectively,  (c)  Final  states 
|k+k^>:  negative-ion  states:  two  electrons  transferred  to  the  ed  level. 

*L  P(»)  versus  9,  the  Incident  angle  of  impact:  v1  <  v2  <  v^  <  v^  . 

5.  PT(«0  vs  X  for  various  values  of  the  temperature  T:  (a)  U  -  1 .5  eV  and  (b)  U 

-  M.6  eV.  The  functional  values  are  given  in  the  inset  for  several  points 

to  illustrate  the  numerical  values. 
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IV.  Summary 

A  perturbative  treatment  using  a  Fano-Anderson  effective  interaction  has 
been  presented  for  one-electron  transfer  processes  In  positive-ion  scattering 
from  surface,  i.e.,  ion  neutralization,  and  the  conditions  under  which  the 
perturbative  series  has  fast  convergence  have  been  discussed.  In  particular, 
fast  convergence  is  obtained  when  the  interaction  time  is  short  enough  and  when 
the  ionic  level  is  not  in  resonance  with  any  states  of  the  valence  band.  For 
two-electron  transfer  leading  to  negative-ion  formation,  the  Anderson 
correlation  energy  U  is  included  in  the  analysis,  and  qualitative  agreement  is 
obtained  with  experiments  on  the  conversion  of  H+  (D+)  to  H  (D  )  by  scattering 
from  a  cesiated  W(110)  surface.  Temperature  effects  are  considered,  and  it  is 
seen  that  negative-ion  formation  is  significantly  enhanced  at  finite  temperature 
T,  provided  kgT  is  not  less  than  U.  Improvements  on  the  calculations  are 

possible  by  incorporating  the  image  potential  into  the  ionic  energy  level  shift. 

9  34*  35 

Finally,  a  more  microscopic  approach  ’  ’is  possible  by  considering  the  ion- 
surface  interaction  strength  V  (see  Eq.  (13)).  Once  this  is  completed,  the 
ion-surface  scattering  can  be  understood  from  a  microscopic  point  of  view. 
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where  eQ  is  the  ionic  energy  when  the  ion  is  far  away  from  the  surface.  The 
decrease  of  the  ionic  energy  as  it  approaches  the  surface  leads  to  the  situation 
that  the  ionic  level  is  in  resonance  with  conduction  electrons  of  the  metal.  As 
such,  electron  tunnelling  is  possible. 

To  make  our  arguments  more  specific,  we  look  at  Eq.  (20),  which  can  be 
written,  with  the  help  of  Eq.  (21),  as 

)  iH'f  -iH’t’f^  iH-t"  -iH-t" 

TUJ(t)  -  (-i r  df  e  u  Hj( t' )  e  u  dt"  e  u  HI(t")e  u  .  (27) 

—  00  —  QD 

Here  we  distinguish  and  Hg,  where  the  former  corresponds  to  the  Hamiltonian 
in  which  one  of  the  electrons  is  already  transferred  to  the  ion,  while  the 
latter  corresponds  to  that  of  the  incident  positive  ion.  They  are  given, 
respectively,  by 


H"(t)  -  l  e.  c*  c,  +  l  e'(t)c*  c.  +  U  n  n. 

0  ,*■  k  ko  ko  *:  d  do  do  do  d,-c 


H'  (t)  -  l  e.c.  c.  +  J  c”(t)c^  c.  +  U  n .  n  (29) 

0  ,  k  ko  ko  1  .  d  do  do  do  d-o 

ko  d 

where  c^(t)  and  e^(t)  are  of  the  form  of  Eq.  (26)  but  with  different  image 
potentials.  The  one  in  e^(t)  is  the  atom-surface  image  potential,  and  that  in 
ejj(t)  is  the  ion-surface  image  potential.  The  forms  of  Hg(t)  and  H^(t)  only 
differ  in  the  ionic  energy  level  due  to  the  change  of  state  resulting  from 
electron  transfer.  Since  the  effect  of  the  image  potential  is  small,  e^(t)  might 
be  taken  as  independent  of  time,  i.e.,  independent  of  position  within  the  framework 
of  the  classical  trajectory  approximation.  We  expect  this  consideration  to  give 
rise  to  better  quantitative  comparison  with  experimental  measurements.^ 
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where  is  the  upper  edge  of  the  conduction  band  and  is  set  equal  to  5  eV  to 

mimic  the  usual  experimental  situation  of  light  atoms  impacting  on  a  transition 
31 

metal  surface.  Performing  calculations  similar  to  the  case  of  T-0,  we  obtain 
the  dependence  of  P^.(»)  on  the  dynamical  variable  A  and  the  temperature  T  for 
different  values  of  U  as  shown  in  Fig.  5.  One  can  readily  see  that  temperature 
effects  are  more  prominent  for  small  U  [Fig.  5(a)],  while  there  is  hardly  any 
effect  for  large  U  [Fig.  5(b)].  One  can  thus  conclude  that  the  negative-ion 
formation  is  significantly  enhanced  at  finite  temperture  T,  provided  knT  is  not 

D 

less  than  the  Anderson  correlation  energy  U. 

In  the  above  calculation,  we  have  assumed  that  the  ionic  level  e  .  is 

d 

constant  in  time.  However,  it  is  well  known  that  a  charge  outside  the  metal 

surface  induces  an  image  of  opposite  charge  inside  the  metal  due  to  the  mobility 

of  conduction  electrons  of  the  metal.  The  induced  image  in  turn  interacts  with 

32 

the  ion  by  an  effective  Coulombic  attraction. 


W 


Mz  ♦  Atf)  ’ 


(2H) 


where  q  is  the  charge  of  the  ion,  z  is  the  distance  between  the  ion  and  the 
surface,  and  is  the  Thomas-Fermi  screening  length  which  is  given  by 

Ayp  *  ^ire2p(ep)  .  (25) 

33 

The  expression  of  Eq.  ( 2H )  has  proved  to  be  a  good  approximation  down  to  the 
distance  of  ATp/2  from  the  surface.  Taking  the  image  potential  of  Eq.  (24)  into 
account,  the  energy  level  for  the  incident  ion  is 

ed(t)  “  G  0+  W 


(26) 
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HjU)  •  exp(iH0t)  Hj  ( t)  exp(-lH0t)  . 


Cl) 


The  experimental  observation  of  outgoing  negative  ions  after  the  scattering 

<=vent  is  complete  is  represented  by  PC®).  Wc  take  =  -10  eV,  f  ^  *  -O.Of  eV,  > 

-  v  cosP  and  A  «  p  \!  -  0 .  ■M  3  exp(-0.01  v)  (in  the  unit  of  eV),  where  v  is  the 

magnitude  of  the  velocity,  P  is  the  angle  between  the  incident  particle  anc  the 

normal  direction  normal  to  the  surface,  and  the  exponential  dependence  of  A  on  v 

accounts  for  the  loss  of  particles  due  to  penetration  into  the  surface.  Here  we 

also  assume  that  the  ionic  energy  level  is  fixed  during  the  scattering  event, 

and  therefore  v  put  c  as  the  zero  energy  level  in  the  calculation.  The 

c 

results  of  the  calculation  of  P  ( “■ )  versus  grazing  angle  P  for  various  incident 
velocities  of  the  ion  are  shown  in  Fig.  it.  They  are  in  qualitative  agreement 

15  4  4-- 

with  experimental  measurements  on  the  conversion  of  H  (D  )  to  H  (D  )  by 

scattering  from  a  cesiated  K(110)  surface.  Further  improvement  of  the 
theoretical  calculations  will  be  discussed  later. 

We  have  also  extended  the  above  unified  point  of  view  to  consider  the  effect 

30 

of  temperature  on  the  negative-ion  formation  at  a  surface.  At  finite 
temperature  T«0,  the  distribution  of  the  conduction  electrons  in  the  solid  is 
governed  by  the  Fermi-Dirac  distribution  given  as 


f(c) 


_  1 _ _ 

1  ~  expiU-tjJ/kgTj  ’ 


(22) 


where  kD  is  the  Boltzmann  constant.  Taking  this  into  account,  the  probability 

D 

of  negative-ion  formation,  P^.(t),  at  time  t  for  the  case  of  finite  temperature 
is 


PT(t) 


i " 


de  p(t)fU)  I  de’  pU’)f(e’)  |  <k+k; |f (2) ( t) 1 1> |2  , 

•  c 


(23) 
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